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Abstract 

Some simple examples from quantum physics and control theory 
are used to illustrate the application of the theory of Lie systems. 

We will show, in particular, that for certain physical models both 
of the corresponding classical and quantum problems can be treated 
in a similar way, may be up to the replacement of the involved Lie 
group by a central extension of it. 

The geometric techniques developed for dealing with Lie systems 
are also used in problems of control theory. Specifically, we will study 
some examples of control systems on Lie groups and homogeneous 
spaces. 

MSC: Primary 34A26; Secondary 81Q70, 93B29. 



1. Introduction: Lie systems. 

There exists a class of systems of time-dependent first order differential 
equations 

(11* 

— = X\x\...,x n ,t), i = l...,n, (1) 

for which there is a function $ : M n ( m + 1 ) — » W n such that the gen- 
eral solution can be written as x = $(xn\, . . . , xi m y, k\, . . . , k n ), where 
i x (j) I 3 = 1> • • • > m } i s an y se t °f particular but independent solutions 
of (pQ) , and ki , . . . , k n , are n constants characterizing each particular 
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solution. These systems, to be called Lie systems, have been charac- 
terized by Lie ^S] an d are receiving much attention in recent years, 
both in physics and in mathematics From the geometric view- 

point, Lie systems correspond to i-dependent vector fields which can 
be written as a linear combination, with i-dependent coefficients, of a 
finite set of true vector fields closing on a finite-dimensional real Lie 
algebra 0[7j. 

A simple instance of Lie system is the linear system 

^ = X>*i(t)^\ i = l,-..,n, (2) 

for which the general solution can be written as a linear combina- 
tion of n independent particular solutions xm, . . . , xr n \, i.e., x = 
3>(cc(i), . . . , xr n \, fci, . . . , k n ) = ki + • • • + k n X( n )i an d in a similar 
way, the general solution for an inhomogeneous linear system can be 
written as an affine function of n + 1 independent particular solutions. 
Another very remarkable example is the Riccati equation 

= a 2 (t)x 2 (t)+a 1 (t)x(t) + a {t) , 

for which the superposition formula comes from the fact that the cross 
ratio of four different solutions is a constant, see, e.g., [TU] . 

The main point is that Lie systems are always related with Lie 
systems on Lie groups defined by right-invariant vector fields. Let G 
be a Lie group. If {ax, ■ ■ ■ , a r } is a basis of the tangent space T e G at 
the neutral element e E G and denotes the right-invariant vector 
field in G such that X^(e) = a a , a Lie system on G will be written as 

g(t) = -j^b a (t)X*(9(t)) ■ (3) 
a=l 

When applying (R g (t)- 1 )*g(t) to both sides we obtain 

r 

(Rg(t)-i)*g(t)(g(t)) = Y1 & o(*) a ° • ( 4 ) 

a=l 

which is usually written, with a slight abuse of notation, as 



a=l 
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This equation is right-invariant: if g(t) is a solution with initial con- 
dition g(0) = e, the solution with initial condition g(0) = go is given 
by g(t)go- Therefore, we only need to find the solution of @ starting 
from the neutral element. 

Let H be any closed subgroup of G, M = G/H the corresponding 
homogeneous space, r : G — > G/H the natural projection, and <1> : 
GxM->M the usual left action of G on M. The right-invariant 
vector fields are r-projectable onto the corresponding fundamental 
vector fields — X a = —X aa , i.e., 

n g X*(g) = -X a (gH) , 

and therefore we have a Lie system on M = G/H associated to (J2J): 

r 

= X}M t )- x a(s(t)) ' 

a=l 

where we denote x = gH € M. Then, the solution of this system 
starting from xo is given by x(t) = $(g(t),xo), where g(t) is the so- 
lution of © starting from the identity. In this sense the equation J3J) 
has a universal character, and it will have an associated Lie system 
on each homogeneous space of G 0. 

Lie systems are of interest not only in the theory of differential 
equations but also in other related fields. For example, they are im- 
portant in classical or even in quantum physics (for instance, in order 
to study the non-relativistic dynamics of a spin 1/2 particle, when 
only the spinorial part is considered |S] ) • Another field where Lie 
systems play an important role is in geometric control theory. 

The aim of this article is to illustrate these applications. Thus, 
after a brief account of a generalization of the method proposed by 
Wei and Norman d HS1 HU, to be used later, we will study 
the particular case where the Lie systems of interest are Hamiltonian 
systems as well, both in the classical and quantum frameworks. The 
theory is illustrated through the particularly interesting example of 
generic classical and quantum quadratic time-dependent Hamiltoni- 
ans. In particular, we show that there exist i-dependent quantum 
systems for which one is able to write in an explicit way the time evo- 
lution of any state of the system. The very simple case of both the 
classical and quantum time-dependent linear potential will be explic- 
itly solved. We will show as well the use of the theory of Lie systems 
in geometric control theory, when dealing with drift-free systems that 
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are linear in the control functions. In particular, we will study from 
this new perspective several well-known control systems: the robot 
unicycle [21,, the Brockett nonholonomic integrator [I] and its real- 
ization in the model of a hopping robot in flight phase j2U|, and the 
kinematic equations of a generalization due to Jurdjevic ^3] of the 
elastic problem of Euler. 



2. The Wei and Norman method. 

Let G be a Lie group as in the previous section. We are interested in 
finding the curve g(t) € G such that 

r 

elegit)- 1 = -Y J b a (t)a a , (5) 

a=l 

with g(0) = e £ G. We can use a method which is a generalization 
of the method proposed by Wei and Norman, in order to find the 
time-evolution operator for linear systems of type 

with U(0) = I. The generalized Wei-Norman method consists of 
writing the previous g(t) in terms of a set of second kind canonical 
coordinates, 

r 

g(t) = exp(-v a (t)a a ) = exp(-«i(t)ai) • • ■ exp(-v r (t)a r ) , 

Q=l 

and transforming the equation © into a system of differential equa- 
tions for the v a (t), with initial conditions v a (0) = 0, a = 1, . . . , r. 
Such a system is obtained from the following relation: 



I TT exp(-u /3 (t) ad(a /3 )) \ a a = ^ b a (t)a a . (6) 

a=l \/3<a ) a=l 

If the Lie algebra of G is solvable, the solution of the previous system 
can be obtained by quadratures. If instead, the Lie algebra of G 
is semi-simple, then the integrability by quadratures is not assured 

n Eid nam. 
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3. Hamiltonian systems of Lie type. 



An interesting and important case occurs when (M, Q) is a symplectic 
manifold and the vector fields in M arising in the expression of the t- 
dependent vector field describing a Lie system are Hamiltonian vector 
fields closing on a finite-dimensional real Lie algebra Q. These vector 
fields correspond to a symplectic action of a Lie group G with Lie 
algebra Q on the symplectic manifold (M, fi). 

The Hamiltonian functions h a of such vector fields, defined by 
i(X a )U = -dh a , in general do not close on the same Lie algebra Q 
when the Poisson bracket is considered, since we can only assure that 

d({h a ,hp} - fya,/?]) = , 

and therefore, they span a Lie algebra extension of the original one. 

The situation in quantum mechanics is quite similar: the Hilbert 
space TL can be seen as a real manifold with a global chart. The 
tangent space TmTL at any point (j) € TL can be identified with TL 
itself, where the isomorphism which associates ip S TL with the vector 
ip € T^TL is given by: 

WW := (j t f{4> + ^)) , V/ e C°°(H) . 
The Hilbert space TL is endowed with a symplectic 2-form defined 

by 

A vector field is just a map A: TL — > TL; therefore a linear operator 
A on TL is a special kind of vector field. Given a smooth function 
a: TL — > M, its differential da^ at € is an element of the (real) 
dual TL' given by 

Actually, the skew-self- adjoint linear operators —i A in TL, for a self- 
adjoint operator A, define Hamiltonian vector fields, the Hamiltonian 
function of —i A being a{(f>) = A<j)). Therefore, the Schrodinger 
equation plays the role of Hamilton equations, because it determines 
the integral curves of the vector field —i H, where H is the Hamiltonian 
of the system . 
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In particular, the theory of Lie systems applies in the previous 
framework when we have a t-dependent quantum Hamiltonian that 
can be written as a linear combination, with t-dependent coefficients, 
of Hamiltonians Hi closing on a finite-dimensional real Lie algebra 
under the commutator bracket. However, note that this Lie algebra 
does not necessarily coincide with that of the corresponding classical 
problem, but it may be a Lie algebra extension of the latter. 



4. Time-dependent quadratic Hamilto- 
nians. 

For the illustration of the classical and quantum situations described 
in the previous section, we consider now the important examples pro- 
vided by the time-dependent classical and quantum quadratic Hamil- 
tonians. 

The first one is the mechanical system for which the configuration 
space is the real line R, the corresponding phase space T*R, endowed 
with its canonical symplectic structure oj = dq A dp, and the time- 
dependent classical Hamiltonian 

H = a(t) £ + (3{t) SI + 7 (t) t + 5(t) p + e(t) q . (7) 

The dynamical vector field solution of the dynamical equation 

i(T H )uj = dH , 

is given by 

(8) 

which can be rewritten as 

= a(t) X 1 + 0(t) X 2 + 7 (t) X 3 - 5(t) X A + e(t) X 5 , 

with 

d 1 / d d \ d 

Xi=p—, X 2 = - U— -p— , X 3 = -q — , 
oq 2 \ oq op J op 



9 X - 8 
oq op 
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being vector fields which satisfy the following commutation relations: 



[X\ , X2) — X\ , \X\ , X 3 ] — 2X2 , 
[X 1 ,X<] = 0, [X l ,X 5 ] = -X 4 , 

[X2, X3] = X 3 , [X2, X4] = — - X4 

[X2,X 5 ] = - X 5 , [X 3 , X4] = X 5 , 
[X 3 ,X 5 ] = 0, [X 4 ,X 5 }=0, 



(9) 



and therefore they close on a five-dimensional real Lie algebra. Con- 
sider the abstract, five-dimensional, Lie algebra g such that in a basis 
{ai, 02, a 3 , 04, as}, the Lie products are analogous to that of ©. 
Then, g is a semi-direct sum of the Abelian two-dimensional Lie alge- 
bra generated by {04, 05} with the sl(2, M) Lie algebra generated by 
{ai, 02, 03}, i.e., Q = M 2 x $1(2, R). The corresponding Lie group will 
be the semi-direct product G = T2 SL(2,M) relative to the linear 
action of SL(2, R) on the two-dimensional translation algebra. When 
computing the flows of the previous vector fields X a , we see that they 
correspond to the affine action of G on R 2 , and therefore, the vector 
fields X a can be regarded as fundamental fields with respect to that 
action, associated to the previous basis of the Lie algebra. 

In order to find the time-evolution provided by the Hamiltonian 
0, i.e., the integral curves of the time-dependent vector field (jSJ, 
we can solve first the corresponding equation in the Lie group G and 
then use the affine action of G on R 2 . We focus on the first of these 
questions: we should find the curve g(t) in G such that 



with 61 (t) = a(t), b 2 (t) = 0(t), 63 (t) = j(t), 6 4 (i) = -S(t), and 
65 (t) = e(t). The explicit calculation can be carried out by using the 
generalized Wei-Norman method, i.e., writing g(t) in terms of a set 
of second class canonical coordinates, for instance, 

g(t) = exp(-w 4 (i)a4)exp(-U5(t)a 5 )exp(-t>i(t)ai)- 
• ex.p(-V2(t)a 2 ) exp(-v 3 (t)a 3 ) , 

and then, a straightforward application of ^ leads to the system 

v 1 = b 1 + b 2 v 1 + 63 v\ , v 2 = b 2 + 2 63 vi , v 3 = e" 2 b 3 , 
1 1 

V4 = 64 + 77 b 2 Vi + 61 V 5 , V 5 = 65 - 63 V4 - - b 2 v 5 , 
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2 



2 
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with initial conditions i>i(0) = ■ • • = i>5(0) = 0. 

For some specific choices of the functions a(t), . . . , e(t), the prob- 
lem becomes simpler and it may be enough to consider a subgroup, 
instead of the whole Lie group G, to deal with the arising system. For 
instance, consider the classical Hamiltonian 

H = ^ + f(t)q, 

which in the notation of (JJJ) has the only non-vanishing coefficients 
a(t) = 1/m and e(t) = f(t). Then, the problem is reduced to one in a 
three-dimensional subalgebra, generated by {Xi, X4., X$}. The asso- 
ciated Lie group will be the subgroup of G generated by {01,04,05}. 
This example will be used later for illustrating the theory: since such 
a subgroup is solvable, the problem can be integrated by quadratures. 

Another remarkable property is that the Hamiltonian functions h a 
corresponding to the Hamiltonian vector fields Xi, . . . , X5, defined by 
i{X a )uj = —dh a , i.e., 

hi(q,p) = h 2 {q,p) = -\qp , h 3 (q,p) = , 

hi(q,p) = P, h 5 (q,p) = -q, 

have almost the same Poisson bracket relations as the vector fields 
X a , but they do not coincide because of {04,05} = 1, instead of 
[X^Xs] = 0. In other words, they close on a Lie algebra which is a 
central extension of M. 2 x st(2, K) by a one-dimensional algebra. 

Let us now consider the quantum case [2^], with applications in 
a number of physical problems, as for instance, the quantum motion 
of charged particles subject to time-dependent electromagnetic fields 
(see, e.g., an d connects with the theory of exact invariants de- 

veloped by Lewis and Riesenfeld (see |18j and references therein). 

A generic time-dependent quadratic quantum Hamiltonian is given 

by 

H = a(t) ^+P(t) QP+ A PQ +lit) ^-+S(t)P+e(t) Q+4>{t)I . (10) 

where Q and P are the position and momentum operators satisfying 
the commutation relation 

[Q,P]=iI . 
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The previous Hamiltonian can be written as a sum with t-dependent 
coefficients 



H = a(t) ffi + (3(t) ff 2 + 7 (i) # 3 - 5(t) ff 4 + e(t) ff 5 - <f>(t)H 6 

of the Hamiltonians 

P 2 1 O 2 

Hi = — , h 2 = -{qp + pq), #3 = ^, 

ff 4 = — P , H§ = Q , Hq = —I , 
which satisfy the commutation relations 



[iff 1 , iff 2 ] = iffi , [iff 1 , iff 3 ] = 2 i ff 2 , [iffi , iff 5 ] = -iff 4 • 

[iff 2 , iff 3 ] = iff 3 , [iff 2 , iff 4 ] = ~ ff 4 , Nff 2 , iffs] = |^5, 
[iff 3 , iff 4 ] = iff 5 , [iff 4 , iff 5 ] = iff 6 



and 

[iffi, iff 4 ] = [iff 3 , 1H5] = [iH a ,iH 6 ] = 0, a = 1, . . . , 5 . 

That is, the skew-self-adjoint operators iH a generate a six-dimensional 
real Lie algebra which is a central extension of the Lie algebra arising 
in the classical case, M? x st(2, R), by a one-dimensional Lie algebra. 
It can be identified as the semi-direct sum of the Heisenberg-Weyl Lie 
algebra f)(3), which is an ideal in the total Lie algebra, with the Lie 
subalgebra sl(2, R), i.e., h(3) xst(2, R). Sometimes this Lie algebra is 
referred to as the extended symplectic Lie algebra f)Sp(2,R) = h(3) x 
Sp(2, R). The corresponding Lie group is the semi-direct product 
ff(3) SX(2,R) of the Heisenberg-Weyl group ff(3) with SX(2,R), 
see also |2*5] . 

The time-evolution of a quantum system can be described in terms 
of the evolution operator U(t) which satisfies the Schrodinger equation 
(see, e.g., PHI) 

i^ = H(t)U, 17(0) = Id, 

where ff (t) is the Hamiltonian of the system. In our current case, the 
Hamiltonian is given by (|10|). and therefore the time-evolution of the 
system is given by an equation of the type 

6 

gg~ l = -Y, h <*( t ) a <* ' 9(0) = e, (11) 

a=l 



9 



with the identification of g(t) with U(t), e with Id, iH a with a a for 
a € {1, . . . , 6} and the time-dependent coefficients b a (t) are given by 

h(t) = a(t), b 2 (t)=(3(t), b 3 (t)=j{t), 
b 4 (t) = -5{t) , 6 5 (t) = e(f) , &e(t) = -0(t) . 

We would like to remark that time-dependent quantum Hamilto- 
nians are seldom studied, because it is generally difficult to find their 
time evolution. However, in the case the system could be treated as 
a Lie system in a certain Lie group, the calculation of the evolution 
operator is reduced to the problem of integrating the system appear- 
ing after the application of the Wei-Norman method. In the case 
the associated Lie group is solvable, the integration can be made by 
quadratures, leading to an exact solution of the problem. We will see 
an example in the next section. 

The calculation of the solution of can be carried out by using 
the generalized Wei-Norman method, i.e., writing g(t) in terms of a 
set of second class canonical coordinates. We take, for instance, the 
factorization 

g(t) = exp(-i;4(t)a4) exp(-u 5 (i)a 5 ) exp(-u 6 (t)o 6 ) 

x exp(-«i(t)ai) exp(-i> 2 (i)a 2 ) exp(-u 3 (t)a 3 ) , 

and therefore, the equation © leads in this case to the system 
i>i = h + b 2 vi + 6 3 v\ , i>2 = b 2 + 2 63 v% , 
V4 = b± + - 6 2 V4 + h v 5 , v 5 = b 5 - 63 Vi - 

ve = be + &5 Vi - - 63 v\ + - bi v\ , 

with initial conditions fi(0) = • • • = v§(f$) = 0. 

Analogously to what happened in the classical case, special choices 
of the time-dependent coefficient functions a(t), . . . , <p(t) may lead to 
problems for which the associated Lie algebra is a subalgebra of that 
of the complete system, and similarly for the Lie groups involved. For 
example, we could consider as well the quantum Hamiltonian linear 
in the positions 

H = £- + f(t)Q , 

which in the notation of (jlUj) has the only non-vanishing coefficients 
a(t) = 1/m and e(i) = f(t). This problem can be regarded as a Lie 



v 3 = e v * 63 , 
1, 

2 b 2 v 5 , 
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system associated to the four-dimensional Lie algebra generated by 
{iHi, iH^, His, iH^}, which is also solvable, and hence the problem 
can be solved by quadratures. 

The treatment of this system, as well as that of its classical version, 
according to the theory of Lie systems, is the subject of the next 
section. 



5. An example: classical and quantum 
time-dependent linear potential. 



Let us consider the classical system described by the classical Hamil- 
tonian 



describing, for instance when f(t) = eEo + eE cos cut, the motion of 
a particle of electric charge e and mass m driven by a monochromatic 
electric field. 

We will study in parallel the classical and the quantum problems 
by reduction of both of them to similar equations, and solving them 
by the generalized Wei-Norman method. The only difference between 
the two cases is that the Lie algebra arising in the quantum problem 
is a central extension of that of the classical one. 



The classical Hamilton equations of motion for the Hamiltonian 
([15)1 are 





and the corresponding quantum Hamiltonian 

„ P 2 




(13) 




(14) 



and therefore, the motion is given by 




(15) 
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can be written as a linear combination 



X = -X 1 -f(t)X 2 , 
m 

with X\ = pd/dq, and X2 = d/dp being vector fields closing on a 3- 
dimensional Lie algebra with X 3 = d/dq, isomorphic to the Heisenberg- 
Weyl algebra, namely, 

[X U X 2 ] = -X 3 , [Xx,X 3 ] = , [X 2 ,X 3 ] = . 

The flow of these vector fields is given, respectively, by 

<M*, (?0,Po)) = (<?o +Pot,Po) , 
4>2(t, (qo,Po)) = (qo,Po + t) , 
4>3(.t,(q ,p )) = (qo + t,p ). 

In other words, {X\, X2, X 3 } are fundamental vector fields with re- 
spect to the action of the Heisenberg-Weyl group H(3), realized as 
the Lie group of upper triangular 3x3 matrices, on M 2 given by 







■ 


l-( 







1 


«1 


03 ) 









1 


0.2 




? 








1 I 







Note that X\, X2 and X 3 are Hamiltonian vector fields with re- 
spect to the usual symplectic structure, $7 = dq A dp, meanwhile the 
corresponding Hamiltonian functions h a such that i{X a )£l = —dh a 
are 



p 2 

h\ = — — , h 2 = q, h 3 = -p 



and therefore 



{hi,h 2 } = -h 3 , {h 1 ,h 3 } = 0, {h 2 ,h 3 } = -l. 

Then, the functions {hi, h 2 , h 3 }, jointly with /14 = 1, close on a four- 
dimensional Lie algebra under the Poisson bracket which is a central 
extension of that generated by {X\, X 2 , X 3 }. 

If {ai, a 2 , a 3 } is a basis of the Lie algebra with non- vanishing 
defining relations [01,02] = —a 3 , the corresponding equation in the 
group H (3) to the system (jri]) reads 

gg' 1 = — ai + f{t)a 2 ■ 
m 
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Using the Wei-Norman formula © with 

g = exp(-ii 3 a 3 ) exp(-ii 2 a 2 ) exp(-ui a\) 
we arrive to the system of differential equations 

iii = — , U2 = ~f(t) , U3 - ui u 2 = , 
m 

together with the initial conditions tii(O) = 1*2(0) = ^3(0) = 0, with 
solution 



t 



f f(t')dt' , «8 = ~ fdt'f 

Jo 771 Jo Jo 



Ui = - , U2 = ~ I f(t') dt' , U 3 = -- I dt' I f(t") dt" 

Therefore, the motion is given by 





( 


1 


±_ 






m 


0- 







1 




\ 













)< 


Po 




K 1 / 



- £ /(*') df 
1 

in agreement with ()15|) , We can immediately identify the constants of 
motion 

h = p(t)+ I f(t')dt' , 
Jo 

h = q(t) - - (pit) + f f(t') dt') t + - f dt' f f(t") dt" . 
m\ Jo J m Jo Jo 

As far as the quantum problem is concerned, notice that the quan- 
tum Hamiltonian H q may be written as a sum 

H q = -Hi- f(t) H 2 , 
m 

with 

p2 

Hi = — , H 2 — —Q . 

The skew-self-adjoint operators —iHi and —iH 2 close on a four- 
dimensional Lie algebra with —i H3 = —i P, and —iH^ = i I, isomor- 
phic to the above mentioned central extension of the Heisenberg-Weyl 
Lie algebra, 

[-iffi,-ijy 2 ] = -iff 3 , [-iffi,-iff 3 ]=0, [-iH 2 ,-iH 3 ] = ^H 4 . 
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-iH q U , 17(0) = Id . 



As we have seen in the preceding section, the time-evolution of 
our current system is described by means of the evolution operator U, 
which satisfies 

dU 
~dt 

This equation can be identified as that of a Lie system in a Lie 
group such that its Lie algebra is the one mentioned above. Let 
{a±, a 2 , 03, 04} be a basis of the Lie algebra with non- vanishing defin- 
ing relations [ai, a 2 ] = a 3 and [a 2 , 03] = 04. The equation in the group 
to be considered now is 

gg' 1 = -— ai + f(t)a 2 ■ 
m 

Using g = exp(— 1*4 04) exp(— 113 03) exp(— u 2 a 2 ) exp(— u\ a±), the Wei- 
Norman method provides the following equations: 

«i = ^ . u 2 = -f(t) , 

{l 3 = -^U2, ill = f(t) U 3 + — ul , 

together with the initial conditions «i(0) = 1*2(0) = 1*3(0) = 1*4(0) = 
0, whose solution is 

«i(t) = - , «2(*) = - f f(t') dt' , u 3 (t) = - f dt' f f(t") dt" , 
m Jo m Jo Jo 



and 



1i4 



These functions provide the explicit form of the time-evolution oper- 
ator: 

U(t,0) =exp(-iu 4 (t))exp(iu 3 (t)P)exp(-iu 2 (t)Q)exp(iu 1 (t)P 2 /2) . 

However, in order to find the expression of the wave-function in a 
simple way, it is advantageous to use instead the factorization 

g = exp(-t> 4 a 4 ) exp(-w 2 a 2 ) exp(-t> 3 a 3 ) exp(-i>i a\) . 

In such a case, the Wei-Norman method gives the system 



«i = £ > v 2 = -f{t) , 

1 
2m 



1 1 2 

t '3 = -m«2, 1>4 = - — «2 
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jointly with the initial conditions i>i(0) = 1)2(0) = v^(0) = 1)4(0) = 0. 
The solution is 

vi(t) = -, v 2 (t) = - f dt' f(t'), 
m f Jo 

V3(t) = ~ f dt' f dt"f(t"), 

m Jo Jo 

^ = -^£ dt '(l! dt " fit "> 

Then, applying the evolution operator on the initial wave-function 
<p(p,0), which is assumed to be written in momentum representation, 
we have 

<t>(p,t) = U(t,0)<t>(p,0) 

= exp(-iv 4 (t))exp(-iv 2 (t)Q)exv(w 3 (t)P) exp(i Vl (t)P 2 /2>(p,0) 

= exp(-iv 4 (t)) eM-iv2(t)Q)e i{v3{t)p+Vl(t)p2/2 U(p, 0) 
= exp(-w 4 (t))e i(,;3W(p+1 ' 2( ' ))+ ^ ( ' ){p ^ 2W)2/2 V(p + v 2 (t),0) , 
where the functions Vi(t) are given by the preceding equations. 




6. Applications in control theory. 

Control systems are described by systems of differential equations 

dx l 

— = F(x\u a ), i = l,...,n, a = l,...,m, (16) 

where u a are the so-called control functions or simply controls, which 
are to be determined in such a way that, e.g., the trajectory passes 
through one or two specific points in the configuration space, or maybe 
gives some cost functional a stationary value. 

A control system is said to be controllable if for any given initial 
point p there exists an integral curve of the corresponding vector field 
along 7r : M. n+m — > M ra such that (ir o j)(0) = p, and a value t\ of the 
parameter of the curve 7 such that (ir o j)(ti) = q for any final point 
<?• 

Consider the case of drift-free systems, linear in the control func- 
tions u a (t), for which the time-dependent vector field, whose integral 
curves are the solutions of (|lfijl. is 

X(t,x) = u x (f)Xi(a;) + • • • + u r (t)X r (x) . 
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Lie systems arise when the vector fields X a close on a finite-dimensional 
real Lie algebra. The cases in which the X a are either right-invariant 
vector fields in a certain Lie group G, or vector fields in a homoge- 
neous space of G, can be dealt with according to the theory of Lie 
systems: reducing the problem to solve an equation on the Lie group 
G, of the form 

r 

9(t) = ^2u a (t)X a (g(t)) . 

a=l 

Controllability of control systems on Lie groups has been analyzed 
by Brockett [Sj and Jurdjevic and Sussmann [TJJ. It can be determined 
by studying algebraic properties of the corresponding Lie algebra g. 

Theorem: A drift-free right invariant system on a connected Lie 
group G is controllable if and only if the Lie algebra generated by 
{Xi, . . .,X r } is Q. 

Proof: If f) is the Lie algebra generated by {X±, . . . , X r }, then the 
Lie algebra of the Lie system is not Q but the subalgebra f). The 
orbit of the neutral element e € G is the subgroup H of G with Lie 
subalgebra f). It is then clear that if f) is a proper subalgebra of G, 
the system is not controllable, while it is so when \) = g. 

As examples of application of the theory of Lie systems to specific 
systems treated in control theory, we will study several well-known sys- 
tems: the robot unicycle |171 121j . the Brockett system termed some- 
times as Brockett nonholonomic integrator 4 1, and a system which 
under certain approximation can be reduced to the former, i.e., the 
model of a hopping robot in flight phase |2()j . Afterwards, we will 
study the kinematic equations of a generalization due to Jurdjevic 
|15j of the elastic problem of Euler, and finally we will briefly show 
how the reduction theory of Lie systems can be applied to two of these 
examples. 

6.1. Robot unicycle or model of an automobile 
as a Lie system. 

Our first example corresponds to the robot unicycle (see, e.g., [T7|). 
Essentially, the same control system arises in a very simplified model 
of maneuvering an automobile |21j . 

The configuration space is K 2 x S 11 , with coordinates (xi, x 2 , X3). 
The control system can be written as 

x\ = b 2 (t) sin £3 , x 2 = 6 2 (t)cosx 3 , x 3 = bi(t) , (17) 
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where b\{t) and b 2 if) are the control functions. Its solutions are the 
integral curves of b\ (t) X\(x) + b2(t) X2(x), where 

9 d d 
A 1 = - — , A 2 = sin x 3 h cos x 3 - — . 

0x3 OX\ 0x2 

The Lie bracket of both vector fields, 

d d 

A~ 3 = [X-y, X 2 ] = COSX3 sinx 3 -— , 

oxy 0x2 

is linearly independent from Xi, X 2 . They satisfy 

[X u x 2 ] = x 3 , [X 2 , X 3 ] = , [X u X 3 ] = -x 2 , 

therefore closing on a Lie algebra isomorphic to se(2). This Lie algebra 
has a basis {a±, a 2 , a 3 } for which 

[ai, a 2 ] = a 3 , [a 2 , a 3 ] = , [ai, a 3 ] = -a 2 . 

Writing the solution of the associated problem in SE(2) as 

g(t) = exp(-vi(t)ai) exp(-v 2 (t)a 2 ) exp(-v 3 (t)a 3 ) , 

then the Wei-Norman method leads to the system 

ii\ = b± , «2 = b 2 cos "li ^3 = &2 sin v\ , 

with vi(0) = v 2 (Q) = v 3 (0) = 0. Denoting Bi(t) = / *&i(s)ds, the 
solution is 

ui(t) = Bi(t), v 2 (t)= [ b 2 {s) cos B^sjds, 

Jo 

v 3(t) = / 62(5) sin Bi(s) (is . 
J 

The action of SE{2) on M 2 x S 4 such that Xi, X 2 , X 3 are the asso- 
ciated fundamental vector fields turns out to be 3>((#, a, 6), (xi, x 2 , x 3 )) = 
(xi,x 2 ,x 3 ) with 

x\ = x\— bcosx 3 — asinx 3 , x 2 = x 2 +6sinx 3 — acosx 3 , x 3 = x 3 — 9, 

where (9, a, b) are the second kind canonical coordinates determined 
by the factorization g = exp(6>ai) exp(aa 2 ) exp(6a 3 ). The composition 
law is 

(0, a, b) {9',a',b') = (9 + 9',a' + a cos 9' + b sin 9',b' -a sin 9' + b cos 9') . 
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Then, the general solution of (|17|) is 

-«2 5 -«3), (^10,^20, 2C30))=(^10 + ^3 COS X 30 + V 2 sinx 30 , 

x 2 o - v 3 sin x 30 + v 2 cos x 30 , x 30 + v\) , 

where v\ = v\(t), v 2 = v 2 (t) and v 3 = v 3 (t) are given above. 

In an alternative way, as the vector fields X 2 and X 3 commute, 
there exist coordinates (2/1,2/2,2/3) such that X 2 = d/dy 2 and X 3 = 
d/dy 3 . For instance, 

2/2 = x\ sin x 3 + x 2 cos x 3 , y 3 = x\ cos x 3 - x 2 sin x 3 , 

which can be completed with y\ = x 3 . Then, 

d d a 
x i = « — ^ — ^2^— • 

dyx dy 2 dy 3 

The control system of interest, whose solutions are again the integral 
curves of the time-dependent vector field b\{t) X\ + b 2 (t) X 2l reads 

yi = h(t), y 2 = b 1 (t)v3 + b 2 (t), y 3 = -bi(t)y 2 . (18) 

Now, the expression of the previous action in terms of the coordinates 
(2/1, 2/2, 2/3) is 

$((0, a, b), (2/1, 2/2, 2/3)) = (2/1 - 0, 2/2 cos 9 - y 3 sin 9 - a cos + b sin 6* , 

2/2 sin 6* + 2/3 cos — a sin — 6 cos 9) , 

and hence, the general solution of (fTH)) is $((— t>i, — i^, — 1*3), (2/10, 2/20, 2/3o)), 
i.e., 

2/1 = 2/io + «i , 

2/2 = 2/20 cos 77i + 2/30 sin i>i + v 2 cos t>i + u 3 sin v% , 
2/3 = 2/30 cos ui - 2/20 sin v\ + u 3 cos v\ - v 2 sin , 

where the v^s are those given above. 



6.2. Brockett nonholonomic control system. 

Another interesting example introduced by Brockett, when dealing 
with problems of optimal control and its relation with singular Rie- 
mannian geometry, is related with the three-dimensional Heisenberg- 
Weyl group H(3), which is the lowest-dimensional non-Abelian nilpo- 
tent Lie group. Such a system is very often considered as one of the 
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prototypical examples relating control theory and extremal problems 
in sub-Riemannian geometry, 

It is the control system in M 3 , with coordinates (x, y, z) 

x = h{t), y = b 2 (t), z = b 2 {t)x - bx(t)y , (19) 

where the functions b\(t) and b 2 {t) are regarded as the controls. The 
solutions of this system are the integral curves of the time-dependent 
vector field b x {t) X x + b 2 (t) X 2 , with 

d d d d 

Xl = a fF ' 2 = ~a — ^ X 7T ■ 

ox oz oy oz 

The Lie bracket 

X 3 = [X u X 2 ] = 2^ 

is linearly independent from X\, X 2 , and the set {Xx, X 2 , X$} close 
on the Lie algebra defined by 

[X 1 ,X 2 ]=X 3 , [X 1 ,X 3 ]=0, [X 2 ,X 3 ]=0, (20) 

isomorphic to the Lie algebra f)(3) of the Heisenberg-Weyl group H(3). 

The Lie algebra fj(3) has a basis {a±, a 2 , a 3 } for which the Lie 
products are 

[oi, a 2 ]=a 3 , [oi, o 3 ]=0, [a 2 ,a 3 } = 0. 

A generic Lie system for the particular case of H(3) takes the form 

R g(t)-u g (t)(9{t)) = -h(t)ax - b 2 (t)a 2 - b 3 (t)a 3 , 

and we are now interested in the one with b 3 (t) = for all t, i.e., 

R g (t)-u g (t)(m) = -h(t) ai - b 2 (t)a 2 . (21) 

Writing the solution starting from the identity of H21j) as the product 
of exponentials 

g(t) = exp(-ui(t)ai) exp(-i> 2 (i)a 2 ) exp(-v 3 (t)a 3 ) , 

and applying the Wei-Norman formula © we find the system of dif- 
ferential equations 

vi = h, v 2 = b 2 , v 3 = b 2 V! , 
19 



with initial conditions v\(0) = 1)2(0) = 1*3(0) = 0. The solution can be 
found immediately: 

vi(t) = / bi(s)ds, v 2 (t) = / b 2 (s)ds, 
Jo Jo 

v 3 (t) = I b 2 {s) [ h(r)drds. (22) 
jo Jo 

The preceding solution can be used in order to find the general 
solution of the given system (|19j) . We only need to find a suitable 
parametrization of the Lie group H(3), and the expression of the group 
action with respect to which the original vector fields are the funda- 
mental vector fields. If we take the canonical coordinates of second 
class defined by g = exp(aa±) exp(ba 2 ) exp(cas), when g £ H(3), it 
can be shown that such an action reads 

$ : H(3) xl 3 — ► M 3 
((a, b, c), (x, y, z)) 1 — > (x — a, y — b, z + ay — bx — ab — 2c) , 

the group law being expressed as 

(a, b, c)(a! , b' , c) = (a + a , b + b' , c + c — ba') . 

Then, the general solution of (|19|) is 

$((-fi, -v 2 , -v 3 ),(xo, yo, zo)) = (xo + vi, y + v 2 , z + x v 2 

-y vi - v^ 2 + 2v 3 ) , 

where v\ = v\(t), v 2 = v 2 (t), and v% = V3(t) are given by (|2*2*|) . 



6.3. Hopping robot in flight phase. 

Next we consider another example coming from a physical model: a 
hopping robot in flight phase, which has been studied, e.g., in |2*U] . 
The system consists of a body with an actuated leg that can rotate 
and extend. The coordinates are (t/j, I, 6), describing the body angle, 
leg extension and leg angle of the robot, respectively. The constant mi 
is the mass of the leg, and the mass of the body is taken to be one. The 
interest is focused on the behaviour of the system for small elongation, 
that is, near / = 0. Precisely, the approximation of the system in 
the neighbourhood of / = will lead to a Lie system related to the 
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Heisenberg-Weyl group, and therefore related also to the previous 
example. 

The controls of the system are the leg angle and extension veloci- 
ties. The control system takes the form |2U| 



i> = bx(t), i = b 2 (t), 9 = ~ ,\, 2 h(t), (23) 



mi(l + l) 2 
l + mi(l + l) 



whose solutions are the integral curves of the time-dependent vector 
field bi(t) Yi + b 2 (t) Y 2 , where now 

d mi {l + l) 2 d Y =- 

1 dij> 1 + mi(l + 1) 2 89 ' 2 dl' 

However, the system ()2.':i|) cannot be considered as a Lie system, since 
the iterated Lie brackets 

[Y 2 , [Y 2 , ...[Y^Yt]---}} 

generate at each step vector fields linearly independent from those 
obtained at the previous stage. Notwithstanding, in order to steer the 
original system by sinusoids, it was proposed in [2U| to take the Taylor 
approximation, linear in /, of the system, that is, 

^ = 6i(t), i = b 2 (t), = -(*i + *aO&i(*) s ( 24 ) 
where the constants k\ and k 2 are defined as 

mi 2m t 
ki = — , k 2 



1 + mi (1+mi) 2 
and then the vector fields become 

= X2 = |. (25) 

Now, the new vector field 

X 3 = [X u X 2 ] = k 2 ^- Q 

closes, jointly with X\, X 2 , the Lie algebra (J50J), so that (|2^)l can be 
regarded as a Lie system with associated Lie algebra f)(3). 
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If we use the previously defined canonical coordinates (a, 6, c) of 
second kind for parametrizing the group g £ H(3), the corresponding 
(local) action to our Lie system reads 

$ : H(3) x M — ► M 
((a, b, c), (ip, I, 8)) i — ► (ip — a, I — b, 9 + k 2 (al — c — ah) + ak±) , 

where M is a suitable open set of R 3 . Then, the general solution of 
the system (|24|) can be written, for t small enough, as 

$((-«l, -«2, -«3),(^o, fo, #o)) = (tpo+vi, lo+v 2 , +k 2 (v3-v 1 l Q -v 1 V2)-k 1 v 1 ) , 

where «i = t>2 = v 2 (t), and V3 = ^(t) are given by l|22|l. This 

result can be checked by direct integration. 

6.4. Kinematics of the generalization of the 
elastic problem of Euler as a Lie system. 

Recently, Jurdjevic has generalized the so-called elastic problem of 
Euler to homogeneous spaces of constant curvature embedded in a 
three dimensional Euclidean space, in order to study certain integrable 
Hamiltonian systems from the point of view of optimal control theory 
|15j . We will only deal with the kinematic equations of such systems, 
wich turn out to be Lie systems. 

The system of interest is the control system with configuration 
space M 3 , and coordinates (x\, x 2 , £3), given by 

Xl = -h(t)x 2 - b 2 (t)x 3 , 

X2 = h(t) Xl + h(t)x 3 , (26) 

£3 = e(b 2 (t)xi - b 3 (t)x 2 ) , 

where e = ±1,0. Its solutions are the integral curves of the time- 
dependent vector field bi(t) X\(x) + b 2 (t) X 2 (x) + b 3 (t) X 3 (x), where 



d 


d 


= Xl n 

ox 2 


ox 1 


d 


d 


dx 3 


Xs dx\ 


d 


d 


= X3 « 

ax 2 


ox 3 



(27) 



22 



These vector fields satisfy the commutation relations 

[Xi, X2] = X3 , [X2, X3] = eXi , [X3, X{\ = X2 , 

and hence they generate a Lie algebra isomorphic to Q t of the Lie 
group G e , given by G = SE(2), Gi = 50(3) and G_i = 50(2, 1). 
Therefore, the case e = essentially reduces to the first example 
studied in this section. We take a basis {a±, a 2 , 03} of g t in which the 
Lie products read 

[ai,a 2 ]=a 3 , [a 2 ,a 3 ] = eai, [o 3 , ai] = a 2 , 

and define the signature-dependent trigonometric functions (see, e.g., 

my. 

e = 1 f sin x e = 1 

C £ (x) = <( 1 e = 5 £ (x) = < x e = 

e = — 1 I sinhx e = — 1 




and 



which satisfy 



e{ > C e (x) ' 



C e {x + y) = C e (x)C e (y)-eS e (x)S e (y) : 
S e (x + y) = C e (x)S e (y) + S e (x)C e (y) , 

as well as C^{x) + eS 2 (x) = 1, and 

dTe(x) - i + e r e 2 (x)- 1 



e v ' 2 (x) ' 

Writing the solution of the problem associated to in the group 
G t as the product g(t) = exp(— vi(t)a{) exp(— Vztyao) exp(— ^3^)03), 
and using the Wei-Norman formula l|B|). we obtain the system of dif- 
ferential equations for vi(t), V2(t) and vs(t): 

vi = bi + eT e (v 2 ){b3 cos «i + 6 2 sint>i) , 

^2 = b 2 cos ui — 63 sin i>i , 
63 cos v\ + 62 sin vi 



V3 



C e {v 2 ) 
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with i>i(0) = v 2 (0) = ^3(0) = 0. Other possible reorderings of the 
factorization in exponentials will give rise to similar systems of equa- 
tions. For e = ±1 the group G e is simple and none of the Wei-Norman 
systems can be integrated by quadratures in a general case. For the 
particular case treated by Jurdjevic ^H] we must put (with our nota- 
tion) 61 (t) = l,b 2 (t) = and 63 (t) = k(t). 

6.5. Reduction of Lie systems in control theory. 

Finally, we would like to point out that there exist a technique for 
reducing the problem of solving a given Lie system in a Lie group 
G to solving a similar Lie system but in a subgroup H, provided a 
particular solution of the problem corresponding to the former in an 
associated homogeneous space G/H is known (see, e.g., 0). This 
reduction procedure can be shown to be useful as well in the study of 
the particular kind of drift-free control systems, linear in the control 
functions, which in addition are Lie systems. 

Take for instance, in the simplified model of maneuvering an au- 
tomobile discussed before, the subgroup H = {(0, 0, b)}. In these 
coordinates, r : SE(2) — > SE(2)/H is t(8, a, b) = (6, a). Taking co- 
ordinates (zx, z 2 ) in M = SE(2)/H, we have that the left action of 
SE(2) on M is given by 

$((0, a, b), (zi, z 2 )) = (zi + 9, z 2 + acoszi + bsinzi) . 

The fundamental vector fields with respect to this action are 

d d d 

Xi =-— , X? = - cos zi — , X3 = -sinzi— , 
az\ 0Z2 0Z2 

which satisfy 

[Xi, X2] = Xf , [X2 , X3] = , [Xi, X3] = -X2 , 

and the equations on the homogeneous space M to be solved are 

zi = —bx(i), z 2 = -b 2 {t) coszi . 

Assume we have a curve on SE(2), gi(t), such that the coordinates of 
its projection r(<7i(t)) = (z\(t), z 2 (t)) satisfy the previous equations. 
For example, we take g± (t) = (zi(t), z 2 {t), 0). Then, we can reduce 
the problem to solving an equation on H, which takes the form 

b(t) =b 2 (t) sin Zl (t), 
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which is just a Lie system for the additive group of the real line. 

As a second and last example, consider again the kinematic equa- 
tions of the generalized elastic problem of Euler treated previously. In 
this case, however, it is advantageous to consider instead of the Lie 
group G € , its universal covering G e , in order to perform the reduction. 
We have that G x = SU_(2), G_i = £17(1,1), and G = SE{2). The 
elements of the group G t can be parametrized by four real numbers 
(a, b, c, d) such that a 2 + b 2 + e(c 2 + d 2 ) = 1, the group law 

(a, b, c, d)(a, b' , c, d') = (a", b" ' , c", d") 

being given by 

a" = aa — bb' — e (cc + dd') , b" = ba + ab' — e (dc — cd) , 



c 



" = ca' + db' + ad - bd' , d" = da 1 - cb' + be' + ad' , 



see also |3J To perform the reduction we take the subgroup H 
generated by a%. The projection T:G e —>G t /H is defined by 



r(a, b, c, d) 



ac — bd bc + ad 



a- 



+ 6 2 ' a 2 + b 2 



Taking coordinates (zi, z-i) in M = G e /H, the left action of G e on M 
reads 

*((o, b, c, d ),(z u z 2 ))=(^-, *g\ , 

where 

JVi = (a 2 - b 2 - e (c 2 - d 2 )) Zl - 2(ab + e cd)z 2 + (ac - 6d)(l -e{z\ + 3%)) 
N 2 = 2(ab - e c£)z\ + (a 2 - b 2 + e (c 2 - (i 2 ))z 2 + (ad + 6c) (1 - e (z 2 + zf )) 
D = a 2 + 6 2 - 2e ((ac + bd)z x + (ad - 6c)z 2 ) + e 2 {c 2 + d 2 )(z 2 + z\) . 

The fundamental vector fields with respect to this action are 
h d d 

X l = Zi— , 

OZ\ OZ2 

x? = -I (1 + e «- 2 |))A_ £2l22 |_, 

which satisfy 

[X 1 , X2] = X% , [X2, X 3 ] = eX ± , [Xi, X 3 ] = -X2 , 
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and the equations on the homogeneous space M to be solved are 

zi = h(t)z 2 - -b 2 (t){l + e {z\ - 4)) - b 3 (t) ez x z 2 , 
z 2 = -h(t) Zl - b 2 {t)e Zl z 2 - h 3 (t)(l - e (zf - z 2 2 )) . 

Assume we have a curve gi(t) on G t , such that its projection r(gi (t)) = 
(zi(t), z 2 {t)) satisfy the previous equations. We can take, for example, 

(1, 0, Zl (t), z 2 (t)) 

9i{t) = 



VI + e(*jf(t) + **(*)) 



Then, we can reduce the problem to an equation on H: if its solution 
is of the form h(t) = (cos(t> (t)/2), sin(u(t)/2), 0, 0), then v(t) satisfies 

v(t) = -h(t) + e (b 3 (t) Zl (t) - b 2 (t)z 2 (t)) . 
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